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Abstract 

A finite range interacting particle system on a transitive graph is considered. 
Assuming that the dynamics and the initial measure are invariant, the normalized 
empirical distribution process converges in distribution to a centered diffusion process. 
As an application, a central limit theorem for certain hitting times, interpreted as 
failure times of a coherent system in reliability, is derived. 
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1 Introduction 



Interacting particle systems have attracted a lot of attention because of their versatile 
modelling power (see for instance [?, ?]). However, most available results deal with 
their asymptotic behavior, and relatively few theorems describe their transient regime. 
In particular, central limit theorems for random fields have been available for a long time 
[?, ?, ?, ?, ?, ?, ?], diffusion approximations and invariance principles have an even longer 
history ([?] and references therein), but those functional central limit theorems that de- 
scribe the transient behavior of an interacting particle system are usually much less general 
than their fixed-time counterparts. Existing results (see [?, ?, ?, ?]) require rather strin- 
gent hypotheses: spin flip dynamics on Z, reversibility, exponential ergodicity, stationar- 
ity . . (see Holley and Strook's discussion in the introduction of [?]). The main objective of 
this article is to prove a functional central limit theorem for interacting particle systems, 
under very mild hypotheses, using some new techniques of weakly dependent random fields. 

Our basic reference on interacting particle systems is the textbook by Liggett [?], and 
we shall try to keep our notations as close to his as possible: S denotes the (countable) set 
of sites, W the (finite) set of states, X = IV s the set of configurations, and {r] t , t > 0} an 
interacting particle system, i.e. a Feller process with values in X. If R is a finite subset 
of S, an empirical process is defined by counting how many sites of R are in each possible 
state at time t. This empirical process will be denoted by N R = {N R , t > 0}, and defined 
as follows. 

N* = (N R (w)) weW , N t R (w) = J2UVt(x)) , 

where l w denotes the indicator function of state w. Thus N R is a N^-valued stochastic 
process, which is not Markovian in general. Our goal is to show that, under suitable 
hypotheses, a properly scaled version of N R converges to a Gaussian process as R increases 
to S. The hypotheses will be precised in sections [2] and [3] and the main result (Theorem 
14.11) will be stated and proved in section [H Here is a loose description of our assumptions. 
Dealing with a sum of random variables, two hypotheses can be made for a central limit 
theorem: weak dependence and identical distributions. 

1. Weak dependence: In order to give it a sense, one has to define a distance between 
sites, and therefore a graph structure. We shall first suppose that this (undirected) 
graph structure has bounded degree. We shall assume also finite range interactions: 
the configuration can simultaneously change only on a bounded set of sites, and its 
value at one site can influence transition rates only up to a fixed distance (Definition 
13.21) . Then if / and g are two functions whose dependence on the coordinates de- 
creases exponentially fast with the distance from two distant finite sets R\ and R2, we 
shall prove that the covariance between f(rj s ) and g(Ct) decays exponentially fast in 
the distance between R\ and R 2 (Proposition [33]). The central limit theorem 14. II will 
actually be proved in a much narrower setting, that of group invariant dynamics on a 
transitive graph (Definition I3.4p . However we believe that a covariance inequality for 
general finite range interacting particle systems is of independent interest. Of course 
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the bound of Proposition 13.31 is not uniform in time, without further assumptions. 

2. Identical distributions: In order to ensure that the indicator processes {I w (rj t (x)) , t > 
0} are identically distributed, we shall assume that the set of sites S is endowed 
with a transitive graph structure (see [?] as a general reference), and that both the 
transition rates and the initial distribution are invariant by the automorphism group 
action. This generalizes the notion of translation invariance, usually considered in 
U 1 ([?] p. 36), and can be applied to non-lattice graphs such as trees. Several recent 
articles have shown the interest of studying random processes on graph structures 
more general than Z d lattices: see e.g. [?, ?, ?], and for general references [?, ?]. 

Among the potential applications of our result, we chose to focus on the hitting time of 
a prescribed level by a linear combination of the empirical process. In [?], such hitting 
times were considered in the application context of reliability. Indeed the sites in R can be 
viewed as components of a coherent system and their states as degradation levels. Then 
a linear combination of the empirical process is interpreted as the global degradation of 
the system, and by Theorem 14.11 it is asymptotically distributed as a diffusion process 
if the number of components is large. An upper bound for the degradation level can be 
prescribed: the system is working as soon as the degradation is lower, and fails at the 
hitting time. More precisely, let / : w i— > f(w) be a mapping from W to R. The total 
degradation is the real- valued process D R = {Df , t > 0}, defined by: 

D t = £ /H^»- 

wew 

If a is the prescribed level, the failure time of the system will be defined as the random 
variable 

T a R = inf{t >0, Df>a}. 

Under suitable hypotheses, we shall prove that T R converges weakly to a normal distri- 
bution, thus extending Theorem 1.1 of [?] to systems with dependent components. In 
reliability (see [?] for a general reference), components of a coherent system are usually 
considered as independent. The reason seems to be mathematical convenience rather than 
realistic modelling. Models with dependent components have been proposed in the setting 
of stochastic Petri nets [?, ?]. Observing that a Markovian Petri net can also be interpreted 
as an interacting particle system, we believe that the model studied here is versatile enough 
to be used in practical applications. 

The paper is organized as follows. Some basic facts about interacting particle systems 
are first recalled in section El They are essentially those of sections 1.3 and 1.4 of [?], 
summarized here for sake of completeness, and in order to fix notations. The covariance 
inequality for finite range interactions and local functions will be given in section [3l Our 
main result, Theorem 14.11 will be stated in section 01 Some examples of transitive graphs 
are proposed in section El The application to hitting times and their reliability interpre- 
tation is the object of section [61 In the proof of Theorem 14.11 we need a spatial CLT for 
an interacting particle system at fixed time, i.e. a random field. We thought interesting to 
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state it independently in section [3 Proposition [7j] is in the same vein as the one proved by 
Bolthausen [?] on Z d , but it uses a somewhat different technique. All proofs are postponed 
to section [3 

2 Main notations and assumptions 

In order to fix notations, we briefly recall the basic construction of general interacting 
particle systems, described in sections 1.3 and 1.4 of Liggett's book [?]. 

Let S be a countable set of sites, W a finite set of states, and X = W s the set of 
configurations, endowed with its product topology, that makes it a compact set. One defines 
a Feller process on X by specifying the local transition rates: to a configuration r\ and a 
finite set of sites T is associated a nonnegative measure ct(t], •) on W T . Loosely speaking, 
we want the configuration to change on T after an exponential time with parameter 

ct, v = °t(v,0- 
CeW T 

After that time, the configuration becomes equal to £ on T, with probability ct(t}, C)/ct,»?- 
Let denote the new configuration, which is equal to ( on T, and to r] outside T. The 
infinitesimal generator should be: 

n/fo) = £ E CT(vX)(f(v c )-f(v))- (i) 

TcS c&v T 

For Q to generate a Feller semigroup acting on continuous functions from X into K, some 
hypotheses have to be imposed on the transition rates or(rj, •)■ 

The first condition is that the mapping rj \— > cy(r], •) should be continuous (and thus 
bounded, since X is compact). Let us denote by ct its supremum norm. 

ct = sup Ct. v - 

It is the maximal rate of change of a configuration on T. One essential hypothesis is that 
the maximal rate of change of a configuration at one given site is bounded. 

B = sup ct < oo. (2) 

x£S Tlx 

If / is a continuous function on X, one defines Af(x) as the degree of dependence of / on 
x: 

A f {x) = sup{ 1/(7?) - /(C)| , 77, C e X and V {y) = ({y) \/y^x}. 

Since / is continuous, Af(x) tends to as 1 tends to infinity, and / is said to be smooth 
if Af is summable: 

111/111 = E < °°- 
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It can be proved that if / is smooth, then ftf defined by (pQ) is indeed a continuous function 
on X and moreover: 

\\m <B\\f\\. 

We also need to control the dependence of the transition rates on the configuration at other 
sites. If y £ S is a site, and T C S is a finite set of sites, one defines 

c T (y) =sup{ ||cr(77i, •) - c T (r} 2 , -)\\ tv , m(z) = m( z ) ^ z ^ V }> 
where || • \\ tv is the total variation norm: 

\\c T (vi, ■) ~ct{V2, -)\\tv = 2 Yl \°t(Vi,0 -c T (V2,0l 

If x and y are two sites such that x ^ y, the influence of y on x is defined as: 

l{%,y) = c T (y). 

We will set j(x, x) = for all x. The influences 7(2;, ?/) are assumed to be summable: 

M = sup ^^7(x, < 00. (3) 

ye 5 

Under both hypotheses (EJ) and ((3l), it can be proved that the closure of Q generates a 
Feller semigroup {S t , t > 0} (Theorem 3.9 p. 27 of [?]). A generic process with semigroup 
{St , t > 0} will be denoted by {rj t , t > 0}. Expectations relative to its distribution, 
starting from i] = r\ will be denoted by E^. For each continuous function /, one has: 

S t f(r ] )=E v [f( Vt )}=E[f(r ]t )\ Vo = r ] }. 

Assume now that W is ordered, (say W = {l,...,n}). Let M. denote the class of all 
continuous functions on X which are monotone in the sense that firf) < /(£) whenever 
r] < £. As it was noticed by Liggett (1985) it is essential to take advantage of monotonicity 
in order to prove limit theorems for particle systems. The following theorems discuss a 
number of ideas related to monotonicity. 

Theorem 2.1 (Theorem 2.2 Liggett, (1985)) Suppose^ is a Feller process onX with 
semigroup S(t). The following statement are equivalent : 

(a) f e M implies S(t)f G M, for allt>0 

(b) Hi < /i2 implies fiiS(t) < fj,2S(t) for all t > 0. 

Recall that \i\ < \ii provided that f fdfti < f fdfi 2 for any f e M. 

Definition 2.2 A Feller process is said to be monotone (or attractive) if the equivalent 
conditions of Theorem \2.1\ are satisfied. 
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Theorem 2.3 (Theorem 2.14 Liggett, (1985)) Suppose that S(t) and Q are respec- 
tively the semigroup and the generator of a monotone Feller process on X . Assume 
further that Q is a bounded operator. Then the following two statements are equivalent: 

(a) ttfg > fQg + gttf, for all f , g E M 

(b) fiS(t) has positive correlations whenever \i does. 

Recall that \i has positive correlation if f fgdji > (J fdji) (f gdji) for any f,gE A4. 

The following corollary gives conditions under which the positive correlation property con- 
tinue to hold at later times if it holds initially. 

Corollary 2.4 [Corollary 2.21 Liggett, (1985)] Suppose that the assumptions of Theorem 
\2.3\ are satisfied and that the equivalent conditions of Theorem \2.3\ hold. Let rjt be the 
corresponding process, where the distribution of r q$ has positive correlations. Then for 
ti < ti < ■ ■ ■ < t n the joint distribution of (r] tl , ■ ■ ■ , rj tn ), which is a probability measure on 
X n , has positive correlations. 

3 Covariance inequality 

This section is devoted to the covariance of f(r] s ) and g(rjt) for a finite range interacting 
particle system when the underlying graph structure has bounded degree. Proposition 13.31 
shows that if / and g are mainly located on two finite sets Ri and R 2 , then the covariance 
of / and g decays exponentially in the distance between Ri and R 2 . 

From now on, we assume that the set of sites S is endowed with an undirected graph 
structure, and we denote by d the natural distance on the graph. We will assume not 
only that the graph is locally finite, but also that the degree of each vertex is uniformly 
bounded. 

WxeS, \{yeS, d(x,y) = 1}\ < r , 

where | • | denotes the cardinality of a finite set. Thus the size of the sphere or ball with 
center x and radius n is uniformly bounded in x, and increases at most geometrically in n. 

\{ y eS,d(x,y)=n}\<^—(r-ir and \{y G S , d(x, y) < n}\ < ^— (r-l) n . 

r — 1 r — 2 

Let R be a finite subset of S. We shall use the following upper bounds for the number of 
vertices at distance n, or at most n from R. 

\{x G S , d(x, R) = n}\ <\{y G S , d(x, R) < n}\ < 2\R\e np , (4) 

with p = log(r — 1). 

In the case of an amenable graph (e.g. a lattice on Z d ), the ball sizes have a subexpo- 
nential growth. Therefore, for all e > 0, there exists c such that : 

\{x G S, d(x,R) = n}\ < \{y G S, d(x,R) < n}\ < ce n£ . 



6 



What follows is written in the general case, using (jlj). It applies to the amenable case 
replacing p by e, for any e > 0. 

We are going to deal with smooth functions, depending weakly on coordinates away 
from a fixed finite set R. Indeed, it is not sufficient to consider functions depending only 
on coordinates in R, because if / is such a function, then for any t > 0, S t f may depend 
on all coordinates. 

Definition 3.1 Let f be a function from S into R, and R be a finite subset of S. The 
function f is said to be mainly located on R if there exists two constants a and (3 > p such 
that a > 0, (3 > p and for all x G M: 



Since /3 > p, the sum Af(x) is finite. Therefore a function mainly located on a finite 
set is necessarily smooth. 

The system we are considering will be supposed to have finite range interactions in the 
following sense (cf. Definition 4.17, p. 39 of [?]). 

Definition 3.2 A particle system defined by the rates Qr(?7, ■) is said to have finite range 
interactions if there exists k > such that if d(x,y) > k: 

1. Ct = for all T containing both x and y , 

2. nf(x,y) = 0. 

The first condition imposes that two coordinates cannot simultaneously change if their 
distance is larger than k. The second one says that the influence of a site on the transition 
rates of another site cannot be felt beyond distance k. 

Under these conditions, we prove the following covariance inequality. 

Proposition 3.3 Assume (Hjj and Assume moreover that the process is of finite range. 
Let R\ and R% be two finite subsets of S. Let (3 be a constant such that (3 > p. Let f and 
g be two functions mainly located on R± and R2, in the sense that there exist positive 
constants Kf, K g such that, 



■/3d(x,R) 



(5) 



A f (x) < K f e~ pd{x ' Rl) and A g (x) < n g e 



/3d(x,R 2 ) 



Then for all positive reals s, t 



sup Cov v (f( Vs ),g( Vt )) < CK f K g (\Ri\ A \R 2 \) e D( - t+s) e~^' p)diRl ' R2) , 



(6) 



where 
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Remark. Shashkin [?] obtains a similar inequality for random fields indexed by Z d . 

We now consider a transitive graph, such that the group of automorphism acts transi- 
tively on S (see chapter 3 of [?]). Namely we need that 

• for any x and y in S there exists a in Aut(S), such that a(x) = y. 

• for any x and y in S and any radius n, there exists a in Aut(S), such that a(B(x, n)) = 
B(y,n). 

Any element a of the automorphism group acts on configurations, functions and measures 
on X as follows: 

• configurations: a ■ T](x) = r)(a' 1 (x)), 

• functions: a ■ f(rj) = f(a ■ 77), 

• measures: J f d(a ■ /1) = f(a-f) d/i. 

A probability measure /x on X is invariant through the group action if a • /i — /j, for any 
automorphism a, and we want this to hold for the probability distribution of rj t at all times 
t. It will be the case if the transition rates are also invariant through the group action. 
In order to avoid confusions with invariance in the sense of the semigroup (Definition 1.7, 
p. 10 of [?]), invariance through the action of the automorphism group of the graph will 
be systematically referred to as "group invariance" in the sequel. 

Definition 3.4 Let G be the automorphism group of the graph. The transition rates 
ct(ti, •) are said to be group invariant if for any a G G, 

c a( r)(a-T],a-C) = c T (v,0- 

This definition extends in an obvious way that of translation invariance on Z d -lattices ([?], 
p. 36). 

Remark. Observe that for rates which are both finite range and group invariant, the 
hypotheses (j2J) and ((3|) are trivially satisfied. In that case, it is easy to check that the 
semi-group {St , t > 0} commutes with the automorphism group. Thus if fi is a group 
invariant measure, then so is [iS t for any t (see [?], p. 38). In other terms, if the distribution 
of r] is group invariant, then that of rj t will remain group invariant at all times. 



4 Functional CLT 

Our functional central limit theorem requires that all coordinates of the interacting particle 
system {r) t , t > 0} are identically distributed. 
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Let (B n ) n >i be an increasing sequence of finite subsets of S such that 

S = (J B n , lim l —f = , (7) 



n=l ' 1 



recall that | • | denotes the cardinality and dB„ = {x G B n , 3y G" B n , d(x,y) = 1}. 

Theorem 4.1 Let /i = 5 V be a Dirac measure where 77 G X fulfills 7](x) = 7](y) for any 
x,y G S. Suppose that the transition rates are group invariant. Suppose moreover that the 



process is of finite range, monotone and fulfilling the requirements of Corollary \2.4\ Let 
(B n ) n >i be an increasing sequence of finite subsets of S fulfilling Then the sequence 
of processes 



t>0}, for n = 1,2,... 



converges in D([0,T]) as n tends to infinity, to a centered Gaussian, vector valued process 
(B(t,w)) t >o t wew with covariance function T defined, for w, w' G W , by 

r M (s,t)(w, w') = Cov M (l w (r) s (x) ), Iw'(r]t (x))) . 

Remark. One may wonder wether such results can extend under more general initial 
distributions. The point is that the covariance inequality do not extend simply by inte- 
gration with respect to deterministic configurations. We are thankful to Pr. Penrose for 
stressing our attention on this important restriction. Monotonicity allows to get ride of 
this restriction. 



5 Examples of graphs 

Besides the classical lattice graphs in Z d and their groups of translations, which are con- 
sidered by most authors (see [?, ?, ?]), our setting applies to a broad range of graphs. We 
propose some simple examples of automorphisms on trees, which give rise to a large variety 
of non classical situations. 

The simplest example corresponds to regular trees defined as follows. Consider the 
non-commutative free group S with finite generator set G. Impose that each generator g 
is its own inverse (g 2 = 1). Now consider S as a graph, such that x and y are connected 
if and only if there exists g G G such that x = yg. Note that S is a regular tree of degree 
equal to the cardinality r of G. The size of spheres is exponential: \{y , d(x, y) — n}\ = r n . 
Now consider the group action of S on itself: x ■ y = xy: this action is transitive on S (take 
a = yx). 

From this basic example it is possible to get a large class of graphs by adding relations 
between generators; for example take the tree of degree 4, denote by a, b, c, and d the 



9 



generators, and add the relation ab = c. Then, the corresponding graph is a regular tree 
of degree 4 were nodes are replaced by tetrahedrons. The spheres do not grow at rate 4 n : 
\{y , d(x, y) — n}\ — 4 • 3 n / 2 if n is even and \{y , d(x, y) — n}\ — 6 ■ S^ 1 ^ 2 if n is odd. 




Figure 1: Graph structure of the tree with tetrahedron cells. The graph consists in a 
regular tree of degree 4 (bold lines), where nodes have been replaced by tetrahedrons. 
Automorphisms in this graph correspond to composition of automorphisms exchanging 
couples of branches of the tree (action of generator a for example) and displacements in 
the subjacent regular tree. 



6 CLT for hitting times 

In this section we consider the case where W is ordered, the process is monotone and 
satisfies the assumptions in Theorem 14.11 the initial condition is fixed and / is an increasing 
function from W to M. In the reliability interpretation, f{w) measures a level of degradation 
for a component in state w. The total degradation of the system in state r\ will be measured 
by the sum ^2 xeBn f(r](x)). So we shall focus on the process = {D^ , t > 0}, where 

= Df n is the total degradation of the system at time t on the set R = B n : 

A (n) = £ f(vt(x)). 

It is natural to consider the instants at which reaches a prescribed level of degradation. 
Let k = (k(n)) be a sequence of real numbers. Our main object is the failure time T n , 
defined as: 

T n = M{t > 0, L>J n) > k(n)}. 
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In the particular case where W = {working, failed} (binary components), and / is the 
indicator of a failed component, then simply counts the number of failed components 
at time t, and our system is a so-called "fc-out-of-n" system [?]. 

Let wq be a particular state (in the reliability wo could be the "perfect state" of an 
undergrade component). Let rj be the constant configuration where all components are in 
the perfect state wq, for all x G S. Our process starts from that configuration 77, which 
is obviously group invariant. We shall denote by m(t) (respectively, v(t)) the expectation 
(resp., the variance) of the degradation at time t for one component. 

Var £> (n) 

m(t) = E[f(r) t (x)) 1 770 = v] , v{t) = lim P5T ■ 

n-*oo I t> n I 

These expressions do not depend on x G S, due to group invariance. 

The average degradation D^/\B n \ converges in probability to its expectation m(t). 
We shall assume that m(t) is strictly increasing on the interval [0, r], with < r < +00 
(the degradation starting from the perfect state increases on average). Mathematically, 
one can assume that the states are ranked in increasing order, the perfect state being 
the lowest. This yields a partial order on configurations. If the rates are such that the 
interacting particle system is monotone (see [?]), then the average degradation increases. 
In the reliability interpretation, assuming monotonicity is quite natural: it amounts to 
saying that the rate at which a given component jumps to a more degraded state is higher 
if its surroundings are more degraded. 

We consider a "mean degradation level" a, such that m(0) < a < m{r). Assume the 
threshold k(n) is such that: 

k(n) = a\B n \+o(^/\B^\). 

Theorem 14. II shows that the degradation process should remain at distance 0(y/\B n \) 
from the deterministic function \B n \m. Therefore it is natural to expect that T n is at 
distance 0(1/ y/\B n \) from the instant t a at which m(t) crosses a: 

t a = inf{t, m(t) = a}. 



Theorem 6.1 Under the above hypotheses, 

V\B^\ {T n - t a ) — ^ JV(0, a 2 a ), 

n— >+oo 

with: 

2 = v(ta) 
° a (m'(t„)) 2 ' 
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7 CLT for weakly dependent random fields 



As in section [U we consider a transitive graph Q = (S,E), where S is the set of vertices 

and E C y}, x, y G S, x ^ the set of edges. For a transitive graph, the degree r 

of each vertex is constant (cf. Lemma 1.3.1 in Godsil and Royle [?]). 

For any x in S and for any positive integer n, we denote by B(x, n) the open ball of S 
centered at x, with radius n: 

B(x,n) = {y G S, d(x,y) < n}. 

The cardinality of the ball B(x,n) is constant in x and bounded as follows. 

sup \B(x,n)\< 2r n = 2e np =: K n , (8) 
xes 

where p = ln(max(r, 4) — 1): compare with formula (Tj]). 

Let Y = (Y x ) xe s be a real valued random field. We will measure covariances between 
coordinates of Y on two distant sets R\ and R 2 through Lipschitz functions (see [?]). A 
Lipschitz function is a real valued functions / defined on W 1 for some positive integer n, 
for which 

y • - \f(x)-f(y)\ „ 

Lip / := sup — - — | < oo. 

We will assume the the random field Y satisfies the following covariance inequality: for 
any positive real 5, for any disjoint finite subsets R\ and R 2 of S and for any Lipschitz 
functions / and g defined respectively on M} Rl ^ and IR'^ 2 ', there exists a positive constant 
Cs (not depending on / g, R\ and R 2 ) such that 

|Cov(/(F x , x G Ri),g(y x , x G R 2 )\ < C 5 Lip/Lip^ (1^1 A \R 2 \)exp(-5d(R u R 2 )) . (9) 



For any finite subset R of S, let Z(R) = ^2 xeR Y x . Let (B n ) n< zfq be an increasing 
sequence of finite subsets of S such that \B n \ goes to infinity with n. Our purpose in this 
section is to establish a central limit theorem for Z(B n ), suitably normalized. We suppose 
that (Y x ) x( zs is a weakly dependent random field according to the covariance inequality 

In Proposition 17.11 below we prove that, as in the independent setting, a central limit 
theorem holds as soon as Var Z(B n ) behaves, as n goes to infinity, like \B n \ (cf. Condition 
(jTTj) below). So the purpose of Proposition 17.21 is to study the behavior of Var Z(B n ). We 
prove that the limit (flU holds under two additional conditions. The first one supposes 
that the cardinality of dB n is asymptotically negligible compared to \B n \ (cf. Condition 
([7|) in section [I]); the second condition supposes an invariance by the automorphisms of the 
group Q, of the joint distribution (Y x , Y y ) for any two vertices x and y. More precisely we 
need to have Condition (fTUl ) below, 

Cov(F z , Y y ) = Cov(Y a(x) , Y a(y) ), (10) 
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for any automorphism a of Q. 

In order to prove Proposition 17.11 we shall use some estimations of Bolthausen [?] that 
yield a central limit theorem for stationary random fields on Z d under mixing conditions. 
Recall that the mixing coefficients used there are defined as follows, noting by Ar the 
ex-algebra generated by (Y x , x G R), 

a k ,i(n) = snp{\F(A 1 n A 2 ) - ¥(A 1 )¥(A 2 )\, A> G A Rl} < k, \R 2 \ < I, d(R 1} R 2 ) > n}, 

for n G N and k, I G N U oo, 

p(n) = sup{|Cov(Zi, Z 2 )\, Zi G L 2 (A{ Pi y), \\Zi\\ 2 < 1, d(px, p 2 ) > n}. 

Under suitable decay of (ay(n)) n or of (p(n)) n , Bolthausen [?] proved a central limit 
theorem for stationary random fields on Z d , using an idea of Stein. In our case, instead 
of using those mixing coefficients, we describe the dependence structure of the random 
fields (Y x ) xe s in terms of the gap between two Lipschitz transformations of two disjoint 
blocks (the covariance inequality above). Those manners of describing the dependence 
of random fields are quite different. As one may expect, the techniques of proof will be 
different as well (see section HJ). 

Proposition 7.1 Let Q = (S,E) be a transitive graph. Let (B n ) n ^ be an increasing 
sequence of finite subsets of S such that \B n \ goes to infinity with n. Let (Y x ) xG s be a 
real valued random field, satisfying (GJ). Suppose that, for any x G S , ~EY X = and 
su Pa;es Halloo < oo. Lf, there exists a finite real number a 2 such that 

Y VarZ(B n ) 2 

hm m~\ = a ' ( U ) 

Z(B ) 

then the quantity —^==1^ converges in distribution to a centered normal law with variance 
a 2 . V 

Proposition 7.2 Let Q = (S, E) be a transitive graph. Let (Y x ) x& s be a centered real 
valued random field, with finite variance. Suppose that the conditions and [TW are 
satisfied. Let {B n ) n be a sequence of finite and increasing sets of S fulfilling (EJ). Then 

V |Cov(F , Y z )\<oo and lim — ^-Var Z(B n ) = V Cov(F , Y z ). 

— ' n-+oo Br, L — ' 

z<=S 1 n| z£S 

8 Proofs 

8.1 Proof of Proposition 13.31 

Let T denote the matrix (7(2;, y)) Xtye s, and let it operate on the right on the space of 
summable series £i(S) indexed by the denumerable set S: 



u = {u(x)) xeS ^Tu= (Tu(y)) 



yes, 



13 



with : 

xeS 

(We have followed Liggett's [?] choice of denoting by Tu the product of u by T on the 
right.) Thanks to hypothesis ((3l), this defines a bounded operator of £i(S), with norm 
M. Thus for all t > 0, the exponential of tT, is well defined, and gives another bounded 
operator of £i(S): 

t n Y n u 

Q MK)u = }^— r . 

n=0 

If / is a smooth function, then A/ = (Af(x)) x& s, is an element of £i{S) . Applying exp(tr) 
to Af provides a control on S t f as shows the following proposition (cf. Theorem 3.9 of 

[?])• 

Proposition 8.1 Assume (j2j) and (02). Let f be a smooth function. Then, 

A Stf <exp(tT)A f . (12) 
It follows immediately that if / is a smooth function then S t f is also smooth and: 

III&/HI < e tM 



because the norm of exp(tT) operating on £i(S) is e tM . 

A similar bound for covariances will be our starting point (cf. Proposition 4.4, p. 34 of 
[?])• 

Proposition 8.2 Assume (j2J) and ([3]). Then for any smooth functions f and g and for all 

t > 0, one has, 

\\S t fg - (SJ)(S t g)\\ < ( E CT J r(exp(rr)A / )(y)(exp(rr)A 9 )(z) dr. (13) 

In terms of the process {rj t , t > 0}, the left member of (TT3l is the uniform bound for the 
covariance between f{rj t ) and g(r] t ). 

\\S t fg - (S t f)(S t g)\\ = sup E v [f(r] t )g(r] t )} - E,[/(^)]E,[ 5 (^)] 

nex 

A slight modification of (TT3l) gives a bound on the covariance of f(rj s ) with g(r)t), for 
< s < t. From now on, we shall denote by Cov,, covariances relative to the distribution 
of {i]t , t > 0}, starting at r] = r\: 

Covr,{f(Vs),9(TH)) = ® v [f(vs)g(vt)\ - EM^EMvt)}- 
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Corollary 8.3 Assume (J2J) and ([3]). Let f and g be two smooth functions. Then for alls 
and t such that < s < t, 



sup 



Oov,(/(^),^))|< E ( E C ^J r(exp(rr)A / )(y)(exp(rr)A 5t „J(^)dr. 

(14) 

Proof of Corollary 18. 3L We have, using the semigroup property, 

®n\f{Vs)g{Th)] = E v [f( Vs )E[9(Vt) I Vs]] = E.lfi^S^givs)} = S a {f3 t . s g){ V ). 

Also, 

^vidiVt)} = S t g(v) = S.(St- a g)(Tj)- 
Applying (fT3l ) at time s to / and S t - S g, yields the result. □ 



In order to apply (1141) to functions mainly located on finite sets, we shall need to control 
the effect of exp(tr) on a sequence (Af(x)) satisfying (jSJ. This will be done through the 
following technical lemma. 

Lemma 8.4 Suppose that the process is of finite range. Let R be a finite set of sites. Let 
u = (u(x)) xe s l> e an element of l\(S). If for all x G S, u(x) < ae~^ x,R \ with a > and 
(3 > p, then for all y G S, 

|(exp(tl»(y)| < aexp(2tMe^ +p)k )e-P d{y ' R) . 

This lemma, together with Proposition 18.11 justifies Definition I3.ll Indeed, if / is mainly 
located on R, then by ( fl~2i ) and Lemma [8T4l S t f is also mainly located on R, and the rate 
of exponential decay (3 is the same for both functions. 
Proof of Lemma 18.41 Recall that 

Tu (y) = ^2u(x)^(x,y). 

Observe that if ^{x,y) > 0, then the distance from x to y must be at most k and thus the 
distance from x to R is at least d(y, R) — k. If u(x) < ae~^ x,R ' then: 

Yu{y) < 2ae pk e- my ' R) - k) M = 2ae w+p)k Me- f3d{y ' R) . 

Hence by induction, 

T n u{y) < Q 2 n e^ +p)kn M n e-< id ^ R \ 
The result follows immediately. □ 

Together with (1141) . Lemma l8~4l will be the key ingredient in the proof of our covariance 
inequality. 



15 



End of the proof of Proposition 13.31 Being mainly located on finite sets, the functions 
/ and g are smooth. By (ITlj) . the covariance of f(rj s ) and g(rj t ) is bounded by M(s, t) with: 

M(s,t)= ( E C ^J [\^Mrr)A f )(y)(eMrr)A St _ s9 )(z)dr. 

y,z€S \T3y,z J J ° 

Let us apply Lemma f8T4l to Af and As t _ a9 . 

(exp(fT)A f )(y) < K f exp(rMe ( ' 3+ '' )fc )e- /J *' Rl) = K f e Dr e- pd( - y ' Rl \ (15) 
The last bound, together with (TT2T) . gives 

A St „ s9 (x) < (exp((t - s)T)A g )(x) < K g e D ^e-W x > R2 \ 

Therefore : 

(exp(rr)A St „ s9 )(^) < K g e D ^ +t -^e~^ R2 \ (16) 
Inserting the new bounds (TT5l) and (TTBD into M(s,t), we obtain 

M(s,t)< X! ( E c T ) K/K9 e-«*' fll » +d ( 2 ' fl2 » f S e D ^ +t - s Ur. 

y,zeS \T3y,z J J 

Now if d(y, z) > k and y,z G T, then ct is null by Definition 13.21 Remember moreover 
that by hypothesis ([2]): 

B = sup c r < oo. 

u€S rr-. 

Therefore : 

M(s,t) < Kf>* 9 2D E e-^ fll ) + ^ 2 ». (17) 

y&S d(y,z)<k 

In order to evaluate the last quantity, we have to distinguish two cases. 

• If d{R u R 2 ) < k, then 

^ ^ e - p{d{y ' Rl)+d(z > R2)) < 2e pk ^2e~ l3d{y ' Rl) 

yeS d(y,z)<k yeS 

nGN y&S 
oo 

(p-/3)n 



< 



e 

n=0 



1 _ e -{/9-p) 

4 (p+/3)fc 

< | fl J 4e e -W,ife) 
- 1 1 ii_ e -0»-p) 



1 _ e -(/3-p) 
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• If d(Ri, R 2 ) > k, then we have, noting that d(y, R\) + d(z, R 2 ) > d(R 1 , R 2 ) — d(y, z) and 
that d(y, z) < k, 



E E • 



y&S d(y,z)<k 



Pd{y,Ri) 

d(y,R 1 )<d(R 1 ,R 2 )-k d(y,z)<k d(y,R 1 )>d(R 1 ,R 2 )-k d(y,z)<k 



K ^ R ^ e p(d(R 1 ,R 2 )-k) e pk e -f3(d(R 1 ,R 2 )-k) + ^ R ^ e pk ^ g (p-/3)n 

n>d(Ri,R 2 )-k 



< 41^1 1 + ^— e^^ 1 '^ 



By inserting the latter bound into (jT7L one obtains, 

M{s,t) < CKfKg^le^+^e-V 3 -^*^ , 

with : 

The covariance inequality ([6l) implies that the covariance between two functions essentially 
located on two distant sets decays exponentially with the distance of those two sets, what- 
ever the instants at which it is evaluated. However the upper bound increases exponentially 
fast with s and t. In the case where the process {rj t , t > 0} converges at exponential speed 
to its equilibrium, it is possible to give a bound that increases only in t — s, thus being 
uniform in t for the covariance at a given instant t. 

8.2 Proof of Theorem ED 
8.2.1 Finite dimensional laws 

Let Q = (S,E) be a transitive graph and Aut{Q) be the automorphism group of Q. Let 
/i be a probability measure on X invariant through the automorphism group action. Let 
{Vt)t>o be an interacting particle system fulfilling the requirements of Theorem 14. 1[ Recall 
that {S t , t > 0} denotes the semigroup and nS t the distribution of r/t, if the distribution 
of 7/o is fi. 

Proposition 8.5 Let (B n ) n be an increasing sequence of finite subsets of S fulfilling Q. 
Let assumptions of Theorem \4-l\ hold. Then for any fixed positive real numbers t± < t 2 < 



• • • < tk) the random vector 
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converges in distribution, as n tends to infinity, to a centered Gaussian vector with covari- 
ance matrix (T^ti, tj))x< id < k . 

Proof of Proposition 18.51 We will only study the convergence in distribution of the 

vector 

the general case being similar. For i = 1,2, we denote by a% = (oci(w)) we w two fixed 
vectors of W w \. We have, denoting by • the usual scalar product, 



\ B n\ 



1 

^ E [E ( E aM(UvtM)-^(vtM = «0)j) 

"I x£B n \i=l \w€W J J 



V l-°n| x . eBn 

where (i^^es is the random field defined by 

F * = E ( E ^H0U%(*)) -VMx) = «,))) =: F 1 ( % (x)) + F 2 ( ?7t2 (a;)). (18) 
i=i / 

The purpose is then to prove a central limit theorem for the sum ^2 xeBn Y x . For this, we 
shall study the nature of the dependence of (Y x ) xe s- 

Let i?i and R 2 be two finite and disjoints subsets of S. Let k\ and k 2 be two real valued 
functions defined respectively on M^ 1 ' and M} R ^. Let i^x, K 2 be two real valued functions, 
defined respectively on W Rl and W R2 , by 

l^M) = kj{F x {y{x)) +F 2 { V {x)), x E Rj), j = 1,2. 

Let C be the class of real valued Lipschitz functions / defined on IR n , for some positive 
integer n, for which 

j ■ f \f(x)-f(y)\ 

Lip / := sup — - — | < oo. 

We assume that fci and fc 2 belong to C. Recall that 

Cov^(A;i(y x , x e Ri),k 2 (Y x , x e R 2 )) = Cov, n (K^rfa, r] t2 ), K 2 {r] tl , %)) 

But ^ _^ 

l-^i (^tn%) - K i(v' tl ,Vt 2 )\ <4Lipkx Mw)\ Y \Vti(x) - rf tl (x)\ 
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Denote = 4Lip k\ ^2 weW \ a ii. w )\- Then, the functions 

r] tl — > (Lipfci)Ai(W) ± Kifa^Tfa) 

are increasing. Hence, the functions 

Gf : — > Lip A* ^ (AiCWOT/Jz) + A 2 (W)7fe(x)) ± K x ( Vtl ,r] t2 ) 

are increasing coordinate by coordinate. This also holds for, 

Gf : ^Lipfc 2 ^(^(^^(i) + A 2 (W) Vt2 (x)) ± K 2 (r) tl ,ri t2 ). 

x£R 2 

Under assumptions of Theorem E31 and of its Corollary l2.4l the vector (77^, r] t2 ) has positive 
correlation so that 

CoY v (Gf(r] tl ,r] t2 ),G^(r] tl ,r] t2 )) > 0. 

This gives 

\Cov v {ki{Y x , x E R 1 ),k 2 (Y x , x G R 2 ))\ 

< LipfciLipfe ^ ^ Cov^A^H/)^^) +A 2 (M/)r 7t2 (x) ) A 1 (H/)r 7tl (y) +A 2 (W)r7 t2 (y)). 

From this bilinear formula, we now apply Proposition 13.31 and obtain the following covari- 
ance inequality: for finite subsets Ri and R 2 of S, we have letting 5 = (3 — p, 

|Cov„ (i<Ti (77 tl , 77 42 ) , (?7ti , ?7t 2 ) ) I < CaLip fciLipfc 2 (|i?i| A |# 2 |)exp {-8d{R u R 2 )) , 

where C,$ is a positive constant depending on (3 and not depending on R 1% R 2l fci and fc 2 - 
We then deduce from Proposition 17.11 that i — Sxes ^ converges in distribution to a 

centered normal law as soon as the quantity Var M (^ ;cgB Y x )/\B n \ converges as n tends to 
infinity to a finite number a 2 . This variance converges if the requirements of Proposition 
17.21 are satisfied. For this, we first check the condition of invariance ( 1X01) : 

Cov^(Y x ,Y y ) = Cov^(Y a{x) ,Y a{y) ), 

for any automorphism a of Q and for Y x as defined by ( 1X51 ). We recall that the initial 
distribution is a Dirac distribution on the configuration r\. Then it has positive correlations. 
We have supposed that rj(x) = T](y) for all x,y G S, hence a-/i = fi and the group invariance 
property of the transition rates proves that \i = 8 V fulfills (TX9l) below and then (1X01) will 
hold. Condition ( TX9l ) is true thanks to the following estimations valid for any suitable real 
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valued functions / and g, 



djj,(r}) a ■ S tl (fS t2 - tl g) (77) since a-fi = fi 

d/j,{r])S tl ((a - f)S t2 _ tl (a - g))(r}) since a ■ (SJ) = S s (a ■ f) 
E„((a • /)(%)(a • g)( Vt2 )) = E„(/(a ■ % )p(a • (19) 



Hence Proposition 17.21 applies and gives 

a 2 = ^Cov M (F ,^) 

2 

X;Cov M (lL(ife(0)), 

i,j=l w,w'£W zGS 
2 

».i=i 

where r M (tj,i,) is the covariance matrix as defined in Theorem l4.lt with this we complete 
the proof of Proposition 18.51 

8.2.2 Tightness 

First we establish covariance inequalities for the counting process. Denote g s ,t,w(Viy) = 
^w(i]t(y)) — ^w(Vs(y)) and for any multi-index y = (yi, . . . ,y u ) G S u , for any state vector 
w = (w 1 ,...,w u ) G W u , n yjW = N2=i9t,t,vit(v>yt)- Following ©, for (3 > p, for any 
r-distant finite multi-indices y G S u and z 6 5" , for any times < s < t < T and for any 
state vectors w G W u and w' G PF" 

|Cov„ (n y , W) n ZjW ,) I < AC(u A i;)e 2DT e- (/3 - rfr = c (u A w)e" cr , (20) 
for c = /3 - p and c = Me ^ (1 _ ^) ^ 

Lemma 8.6 There exist 8 > anc? > such that for \s — 1\ < S : 

|Cov„(n X)W ,n yiW 0| <K a \t-s\. (21) 

Proof. Denote f{rf) = 1^(77(2;)) then g t+hAw (r],x) = S h f(r] t ) - f{i] t ); the properties of the 
generator Q imply that 

lim M!tfi = Q/( , ) 



20 



But 

\wm < EE <*{v,o\ftf)-nv)\ 

TdSC,&W T 

< Yl CT ^) < Yl c T<Cn 

TcS,xET TcS,x£T 

so that for h > tending to zero 

\g a ,a+h,v>('ni x )\ - C n h + o(h) 

Because O is group invariant, the remainder term is uniform with respect to index x, so 
that we find convenient 5q and Kq uniformly with respect to location. □ 



From inequality (|20j) and lemma [8761 we deduce the following moment inequality: 

Proposition 8.7 Choose I and c such that p(2l — 1) < c. For (s,t) such that \t — s\ < 
5 Ac e c /K n : 

E{N fn_ N B n) 2i<(U-mcoe 2c )* 



(2/)!(2/-i; 



Cl 

where c\ = p A (c — p(2Z — 1)). 

Proof. Recall that N t B " - iVf™ = Y^xeB n 9s,t,w(v, x). Note that the value of n x 

does not depend on the order of the elements Xi, . . . , xl. The index x is said to split into 
y = (yi, Um) and z = (z ± , z L _ M ) if one can write y 1 = x a{1) , ...,y M = £<x(m) and 
%i = Xcr(M+i), ■ ■ • > zl-m = Xa(L) for some bijection a : {1, . . . , L} — > {1, . . . , L}. We adapt 
lemma 14 in Doukhan & Louhichi [?] to the series (gt, s ,w{Vi x ))xeB n - For any integer q > 1, 
set : 

= £ |En XiW | , (23) 

then, 

E(iVf" - N t Bn ) 21 < \B n \- l A 2l (n). (24) 

If q > 2, for a multi-index x = (x±, . . . ,x q ) of elements of S, the gap is defined by the 
maximum of the integers r such that the index may split into two non-empty sub-indices 
y = (yi, . . . , y h ) and z = . . . , z q -h) whose mutual distance equals r: d(y(x), z(x)) = 
min{d(y a , z;,); l<a</i,l<6<g — /i}=r. If the sequence is constant, its gap is 0. 
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Define the set G r (q, n) = {x G B% and the gap of x is r}. Sorting the sequences of indices 
by their gap: 

n 

A q (n) < Yl E \9s,tAv,xi)\ q + J2 E | Cov ( n y(x),w,n z(x)jW )| (25) 

xieB n r=l xeG r (q,n) 

n 

+ E E |E(n y(x)iW )E(n z(x)iW )|. (26) 

r=l xeG r (q,n) 

Denote 

n 

V q (n) = + E E l Cov ( n y(x) 

n Z ( X ) iW ) I . 

xieB n r=l xeG r (g,n) 

In order to prove that the expression ([26]) is bounded by the product A h {n)A q ^ h {n) we 
make a first summation over the x's such that y(x) e B%. Hence: 

9-1 



A q (n) < V q (n) + J2A h (n)A q . h (n). 



h=l 



To build a multi-index x = (x%, . . . , x q ) belonging to G r (q, n), we first fix one of the \B n \ 
points of B n , say X\. We choose a second point x 2 with d(xi,x 2 ) = r. The third point x 3 
is in one of the ball with radius r centered in one of the previous points, and so on. . . Thus, 
because the maximal cardinality of a ball with radius r writes b{r) < e pr 

\G r (q,n)\ < \B n \b(r)2b(r)---(q-l)b(r) < |5 n |(g - 1)!2^V^. 

We use lemma [8^61 to deduce: 

/ oo 

V q (n) < \B n \ [Knlt-sl + iq-iy^J^^-^icoqe-^AKnlt-sD 



r=l 



Let R be an integer to be specified, then 

(R-l oo 
K u \t - s\ eP(q ~ 1)r + c o E eipi "'~ ' 
r=0 r=R 

Comparing those summations with integrals: 

V g (n) < \B n \q\2^ ( K f~ s L p^-m + c ^_M^)-c)(R-A 

\p(q-l) c-p(g-l) J 

q -i Ka\t - s\ Mq _ 1)R ( c nC _ cR 



< \ Bri \q\2^^ LgPW-^- 1 + -e 
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where c\ = p A (c — p(2l — 1)). Assume that (s, t) G T are such that |£ — s\ < c e c /Kn and 
choose R > 1 as the integer such that e c ^ -1 ^ < „ c °,f < e cR . 

.MJt-s e 2 ^ 1 ' / c 



V» < j , (27) 

so that V q (n) is a function of q that satisfies condition (7i ) of Doukhan & Louhichi [?]. 
Then 

Mn) < _g^_(^(„) + V a („)') 

(4/-2)!(c e 2 ^ / 2 2 '(2Q!(c e 2 -)^ ^ 
" (2/)!(2/-l)! c, \Bn\{K a \t-s\) - 



+ 



Cl 



and Proposition 18.71 is proved. □ 
To prove the tightness of the sequence of processes N Bn , we study its oscillations: 

w (5, N Bn ) = sup | N t Bn - Nf n | 

||t-s||i<<5 

Fix £ and r\. We have to find 5 and n such that for all n > n : 

P(w(5,iV B ") > e) < r] 

Define n as the smallest integer such that \B no \ > <5 _1_p / c , then for n > n , \t — s\ < 5, 
I = 2 and c > 3p, Proposition 18.71 yields: 

E(jV t B ™ - iVf") 4 < C^-c) 

and we now follow the proof in Billingsley [?] to conclude. 



8.3 Proof of Theorem EH 

The proof is close to that of the analogous result in [?]. The convergence in distribution 
of Z n = (Z n (t))t>o, where Z n (t) = (D^ — \B n \ ■ m(t))/ y/\B n \, does not directly imply the 
CLT for T n . The Skorohod-Dudley-Wichura representation theorem is a much stronger 
result (see Pollard [?], section IV. 3). It implies that there exist versions Z* of Z n and 
non-decreasing functions (fi n such that for any fixed s such that for Z*, limit in distribution 
of Z n \ 

lim sup \Z*(t) - Z*((j) n (t))\ = a.s. 
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and: 

lim sup \4> n (t) — t\ = a.s. 

n ^°°0<t<s 

Since Z* has continuous paths, it is uniformly continuous on [0, s], and hence: 

lim sup \Z*(t) -Z*(t)\ = a.s. , (28) 

n^oo 0<t<s 



We shall first use ( 1281) to prove that the distributions of \/\B n \(T n —t a ) are a tight sequence, 
tive constant. On the one hand, if D^ r " 

c/\J\B n \. Thus: 



Let c be a positive constant. On the one hand, if D , > k(n), then T n < t a + 

t a +c/W\B n \ 



nV\B7\{T n - t a ) < c] > ¥[Dl n) > k(n)} 



¥[Z* n (t a + c/v^Bj) > v^J(« - m{t a + c/v^Bj)) + o(l)] 



= FlZ* n (t a + c/^/\B^\) > -cm'(t a ) + o(l)] 
= P[Z*(t Q ) > -cm'(g]+o(l) , 

using (1281) and the continuity of Z*. Since m'(t a ) > 0, we obtain that: 

lim liminf F[x/\B^\(T n - t a ) < c] = 1. (29) 

On the other hand, we have: 



c— >oo n— >oo 



P[v^(T n - t a ) < -c] = ¥[3t <t a - c/ v^Bj , > v^J(« - ™(*)) + o(l)]. 

But since the function m is increasing, for all t < t a — c/ \/\B n \ we have: 



V^J(a - m(t)) > v^Bj(a - m(t a - c/y/\EQ)) = cm'(t a ) + o(l). 

Hence: 



P[v^Bj(T„-t Q ) < -c] < F[3t<t a -c/^/\B^\,Z:(t)>cm'(t a ) + o(l)} 

< P[3i < t„ , Z*(i) > cm'(f Q ) + o(l)] 
= P[3i < * a , Z*(£) > cm'(t a ) + o(l)] + o(l). 

The process Z being a.s. bounded on any compact set and m'{t) being positive on [0,r], 
we deduce that: 

lim lim sup P[ \/\EQ{T n - t a ) < -c] = 0. (30) 



Now ([291 and ( 1301) mean that the sequence of distributions of (\/\B n \(T n — t a )) is tight. 
Hence to conclude it is enough to check the limit. 
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Using again (1281) . together with the almost sure continuity of Z yields: 



D 



' t a +c/y/\R^\ 



\B n \m(t a + u/^\l%\) + y/\B^\Z*(t a + u/^\l%\) + o(^B~T) a.s. 
\B n \a + u^/\^\m\t a ) + ^/\B^Z*(t a ) +o(^/\B^\) a.s. 



Therefore: 



inf { u ; D (n) , > k(n 



inf \u ■ u^\L\\m'(t a ) + y/\B^\Z*(t a ) + o(y/\B^\) > 
Z*(t a 



m'(t a ) 



+ o(l). 



The distribution of —Z*(t a )/m'(t a ) is normal with mean and variance a 2 al hence the 
result. 



8.4 Proof of Proposition 17.11 

Let JF 2 3 be the set of real valued functions h defined on M, three times differentiable, such 
that h(0) = 0, ||/i"||oo < +oo, and ||/i^||oo < +oo. For a function h G ^2,3, we will denote 
by 62 and 63 the supremum norm of its second and third derivatives. We first need the 
following lemma. 

Lemma 8.8 Let h be a fixed function of the set -7-2,3- Let R be a fixed and finite subset 
of S. Let r be a fixed positive real. For any x G R, let V x = B(x, r) n R. Let (Y x ) xe s be 
a real valued random field. Suppose that, for any x G S, EY X = and EY 2 < +00. Let 

z (R) = 12 x eR Y v Then 



E(h(Z(R))) -VarZ(R) / tE(h" (tZ(R)))dt 

Jo 

< [ y2\ c °v(Y xl h\tz(v:)))\dt+2y2nYx\\z(v*)\ib2Ah\z(v x 



^2(Y x z(y x )-E(Y x z(y x ))) 



x€R 



+ b 2 J2\Cov(Y x ,Z(V x c ))\ 



(31) 



xGR 



where V x c = R\V X . 



Remark. For an independent random field (Y x ) xe g, fulfilling sup xe5 EFj 1 < +00, Lemma 
18.81 applied with V x = {x}, ensures 



E(h(Z(R)))-VaxZ(R) \ ffi(h" (tZ(R)))dt 





< 2 E l^! 2 ( fo 2 A &3|n|)+&2Vl^l SUP ||l; 2 ||2 



x£R 



x€S 
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Proof of Lemma 18.81 We have, 



h{Z{R)) = Z{R) [ h'{tZ(R))dt= ! (y2Y x ti(tZ(R)))dt 
Jo Jo \xeR J 

= f (j2Y x h'(tz(v x c )))dt+ f [Y,YAh'{tz{R))-h l {tz{y:))-tz{y x )h''{tz{R))) 

Jo W J Jo \xeR 

+ Y,Y X Z{V X ) [ tti'(tZ(R))dt-J2^( Y x Z (V x )) [ th"(tZ(R))dt 
x&R Jo x & R Jo 

+ ^2E(Y X Z(V X )) I th"(tZ(R))dt-y2E(Y x Z(R)) [ th"(tZ(R))dt 

XER J ° XER J ° 

+ J2E(Y x Z(R)) [ th"{tZ(R))dt. 

xtR J ° 

We take expectation in the last equality. The obtained formula, together with the following 
estimations, proves Lemma [8.81 

\h\tZ{R)) - h'(tZ(V^)) - tZ{V x )h"{tZ{R))\ 

< \h'(tZ(R)) - h'(tZ(V:)) - tZ(V x )h"(tZ(V:))\ + \Z(V x )\\h"(tZ(R)) - h"(tZ{^))\ 
<2\Z(V x )\(b 2 Ab 3 \Z(V x )\). □ 



Our purpose now is to control the right hand side of the bound (I3TT) for a random field 
(Y x ) x( zs fulfilling the covariance inequality ([9j) and the requirements of Proposition 17.11 

Corollary 8.9 Let h be a fixed function of the set T 2 $- Let R be a finite subset of S. For 
any x e R and for any positive real r, let V x = B(x,r) fl R. Let (Y x ) xe s be a real valued 
random field, fulfilling the covariance inequality (0J). Suppose that, for any x G S, EY X = 
and sup x€S ||l^||oo < M, for some positive real M. Recall that Z(R) = ^2 xeR Y x . Then, 
for any 5 > 0, there exists a positive constant C(5 } M) independent of R, such that 



sup 



E(h(Z(R)))-VaiZ(R) / tE(h" (tZ(R)))dt 



1/2 



< C{5, M) I b 2 \R\e- 5r + b 3 \R\K r + b 2 \R\ l/2 n r ^ K k e~^ k - 2r) 



fc=[3r] 



. 1/2' 
[3r]+l x 



-b 2 \R\ 1/2 ^r £ e~ Sk 



K k 



k=l 



recall that sw\) x&s \B(x,n)\ < K n . 
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Proof of Corollary [8791 

We have 

Hence 



V x c = {yeS, d(x,y)>r}nR. 



d({x},V x c )>r. 
The last bound together with (jHJ), proves that 

^\Cov(Y x ,h'(tZ(V x c )))\ < Csb 2 J2(K\ Al)e-™ 

< C s b 2 \R\e- Sr . 

In the same way, we prove that 

b 2 J2\Cov(Y x ,Z(V x c ))\ < C s b 2 \R\e 



-St 



x£R 



Now 



x&R 



J2E\Y x \\Z(V x )\(b 2 Ab 3 \Z(V x )\) < 6 3 M|i?|supE|Z(V;)| 2 

< b 3 M\R\K r sup ^|Cov(y y ,y z ; 



yes 



The last bound is obtained since |T4| < «r and sup^ g5 ^2 ze s |Cov(l^, < oo (the 
of the last inequality is done along the same lines as that of Proposition 17.21 ) . 
It remains to control 



E 



J2(y*Z(V x )-E(Y x Z(V x ))) 



xeR 



For this, we argue as Bolthausen [?]. We have 



E 



J2( Y xZ(V x ) - E(Y X Z(V X ))) 



Va,r(J2Y x Z(V x )) 



x£R 



x£R yeR 



Hence, since V x C B(x,r) 



E 



J2(Y X Z(V X ) - E(Y X Z{V X ))) 



^E E E E icov(y^,r^ 

iefl x'eB(x,r) yeR y'eB(y,r) 



x&R 

We have, 

ICov^y^YyY^)! < |Cov(K c K c /, K^Ky')! ^d(a;,y)>3r + ICov^^', l^fy) | ld(x,y)<3r- 
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We begin by controlling the first term. The covariance inequality (J9j) together with 
elementary estimations, ensures 

oo 

\Cov(Y x Y x /, Y y Yy/) \ ld(x,y)>3r < ' GoV ( Y x Y x> , Y y Y y>) \ I k <d(x,y)<k+1 

k=[3r] 

oo 

< 2M 2 C 5 e- 6 ^' x '^^l k < <Kx>y)<k+1 

k=[3r] 
oo 

< 2M 2 C 5 e-^ k - 2r \^ y)<k+1 , 

k=[3r] 

the last bound is obtained since, for any x' G B(x,r) and y' G B(y,r), we have, 
d({x, x'}, {y, y'}) + 2r> d({x, x'}, {y, y'}) + d(x, x) + d(y, y) > d(x, y). 

Hence, 

E E E E icMn*v,nnoi 

z£-R i'eB(i,r) y£R y'<=B(y,r) 
oo 

<2M 2 C^ 2 £ EE e ^Wi) 

fc=[3r] z£R yG-R 
oo 

fc=[3r] 

We now control the second term in ( |36l ). Inequality ([9)) and the fact that 
d({%}, {%', V, y'}) < d({x}, {x'}), ensure 

\Cov(Y x Y x >, Y y Y yl )\ ld{x,y)<3r 

< \Cov(Y x ,Y x ,Y y Y y ,)\l d{Xty) < 3r + ICov^.,^,)! \Cov(Y y ,Y y ,)\l d{x>y) < 3r 

<2M 2 c 5 e- Sd ^^'>y>y'»i d{x , y) < 3r . 

We deduce, using the last bound, that 

\Cov(Y x Y x >, YyY y ,) I I d ( x ,y)<Zr 
[3r]+l 

< |Cov(y r K c /, Y y Yyi)\ Id(x,y)<3r^k-l<d({x},{x',y,y'})<k 
k=l 

[3r]+l 

< 2M 2 C$ e ^ 1 ' ) ld(i',j/)<3rld({x},{i'',j/,y'})<fc- 

fc=l 

We have 

Id({a;},{a;',y,j/'})<A: < ^d({x},{x'})<k + Irf(M,{s/})<fc + ^({2}, {?/})<&• 
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Hence, we check that, 

E E E E ^d(x,y)<3rh{{x},{x',y,y>})<k < 3\R\nl r K k . 
x€R x'&B(x,r) y£R y'£B(y,r) 

We obtain combining fl38l) and (l39l) . 

E E E E \Cov(Y x Y x ,,Y y Y y ,)\I d{ 

x,y)<3r 

x£R x'eB(x,r) y£R y'&B(y,r) 



(39) 



[3r]+l 



< 



6e s M 2 C s \R\Kl r e~ 6k K k . (40) 



k=i 



We collect the bounds j35]), ([3TJ) and ([401) , we obtain, 

E 



J2(YxZ(V x ) - E(Y X Z(V X ))) 

x£R 

< C(5, M)\R\ 1/2 



1/2 



J2 K fe +ie-^ | +,-,,„ 

. k=[3r] 



[3r]+l 

E e ~ 5 ^ 
fc=i 



1/2' 



• (41) 



Finally, the bounds ((32]), ((33]), {Ml, {H]), together with Lemma EH] prove Corollary El 

□ 



End of the proof of Proposition 17. ll We apply Corollary EH to the real and imaginary 
parts of the function x — > exp(iux/ y/\B n \) — 1. Those functions belong to the set -7-2,3, 



with fc 2 



u 



\Br, 



and 





11 


3 


\B n 


3/2 



We obtain, noting by (j) n the characteristic function of the normalized sum Z(B n )/ >J\B ri 

VarZ(B n ) 2 f 1 



U« - 1 + 



-it / t(f) n (tu)dt 

\£>n\ Jo 



1/2 



< C(5, M, u) I e- 5r + —f^= + 



[3r]+l 



\B n \ 



1/2' 



I -Br, 



, fc=[3r] 



fc=l 



Let S be a fixed positive real such that 5 > 12p, recall that 

sup \B(x,r)\ < 2e rp =: K r . 
xes 
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Hence 
„(u) - 1 



< C{5,M,u) { e~ Sr + 



jrp 



+ 



J\B r .\ 



1/2 



E 



-(S-p)k 



3pr 



< C(M,p,S,u) I e~ Sr + 



3rp 



+ 



Jfc=[3r] 
e -(5-5p)r/2" 



I -Br> 



[3r]+l 

E ■ 

fc=i 



1/2' 



-(5-p)fe 



For a suitable choice of the sequence r (for example we can take r = |ln|£> n |), the right 
hand side of the last bound tends to an n tends to infinity: 



lim 

n— >oo 



t(j) n (tu)dt 



0. 



(42) 



We now need the following lemma. 

Lemma 8.10 iet cr 2 &e a positive real. Let (X n ) be a sequence of real valued random vari- 
ables such that sup ngN EX 2 < +oo. Let <fi n be the characteristic function of X n . Suppose 
that for any h6I, 



lim 



4>n(u) 



Then, for any ii6l, 



lim < 

n— >+oo 



i n [u) = exp( 



t(p n {t)dt 



u 2 a 2 , 



0. 



(43) 



Proof of Lemma I8.10L Lemma [8.101 is a variant of Lemma 2 in Bolthausen [?]. The 
Markov inequality and the condition sup ngN EX 2 < +oo imply that the sequence (/i n )neN 
of the laws of (X n ) is tight. Theorem 25.10 in Billingsley [?] proves the existence of a 
subsequence /i nfe and a probability measure /i such that /i nfc converges weakly to /i as k 
tends to infinity. Let <p be the characteristic function of \i. We deduce from (1431) that, for 
any u G R, 



1 + ^ 



t(p{t)dt = 0, 



or equivalently, for any u G 



fr'(u) + a 2 ud){u) = 0. 



We obtain, integrating the last equation, that for any u G 



[u) = exp( 



a 2 u\ 



The proof of Lemma [8.101 is completed by using Theorem 25.10 in Billingsley [?] and its 
corollary. □ 

Proposition 17.11 follows from (flTi ). ( 1421) and Lemma 18.101 □ 
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8.5 Proof of Proposition 17. 2L 

We deduce from ([9]) that for any positive real 5 there exists a positive constant Cs such 
that for different sites x and y of S, 

\Cov(Y x ,Y y )\<C 5 e- Sd ^. (44) 

Hence, the first conclusion of Proposition 17.21 follows from the bound ( 1441) . together with 
the following elementary calculations, for p < 5, 

J2\Cov(Y ,Y z )\ < C 5 J2 e M-5d(0,z)) 

oo 

< C 5 ^^exp(-5c/(0,z))I r < d (o, 2 )< r+ i 

zeS r=0 

oo 

< C s ^2 exp(-5r) ^ h(o,z)<r+i 

r=0 zeS 
oo 

< C 5 ^exp(-5r)|5(0,r + l)| 

r=0 

oo 

< C(5,p)J2eM-(S-P», (45) 

r=0 

where C(8,p) is a positive constant depending on 5 and p. 

We now prove the second part of Proposition 17.21 Thanks to ([7]) , we can find a sequence 
u = (u n ) of positive real numbers such that 

I I 

lim u n = +oo, lim -. exp(pu n ) = 0. (46) 

n— >+oo n^+oo -D n 

Let (d u B n ) n be the sequence of subsets of S defined by 

d u B n = {s e B n : d(s, dB n ) < u n }. 

The bound (Jlj) gives 

\d u B n \ < 2\dB n \e u " p , 
which together with the suitable choice of the sequence {u n ) ensures 

lim = 0, (47) 

we shall use this fact below without further comments. Let B% = B n \d u B n . We decompose 
the quantity Var S n as in Newman [?]: 

1 r Var S n = t^-| ^ COV Y v) = T ±,n + T 2,n + T 3tn , 

71 
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I B r , I I B r , 



where 



tbiE E cov(y„n) 

\ JD n\ x 



zeB^ y£B n \B(x,u n ) 



T 2 , n - 



^E E Cov ( y -n) 
T ^ = jh E E Cov ( y -^)- 



x€d u B n y£B n 

Control of T lin . We have, since < \B n \ and applying (jH 



\T hn \ < sup ^ |Cov(n,i;)| < C 5 sup exp(Sd(x,y)). (48) 

3:65 y£S\B(a;,u„) xeS yeS\B(x,n) 

For any fixed i 6 S, we argue as for (1451) and we obtain for p < 6, 



ex V {-5d{x, y)) < C{8) £ exp(-(5 - p)r) < C(<f, p) exp(-(5 - pK) (49) 

y€S\B(x,n) r=[u„] 



We obtain, collecting (I48I1 . (1491 together with the first limit in (14611 : 

lim T ln = 0. (50) 

n^+oo ' 

Control of T 3 n . We obtain using 



\T 3n \ < 

1 - LB. 



^sup£|Cov(n,ig|. (51) 



y£S 



The last bound, together with the limit (l47l ) gives 

lim T 3 , n = 0. (52) 



Control of T 2n . We deduce using the following implication, if x G B% and y is not 
belonging to B n then d(x, y) > u n , that 



T2,n = ^E E C ° V ( y -^) 

We claim that, 



xgB" j£B(i,%) 



^ Cov(y;,ig= ^ Cov(r ,y;), (53) 

y£B(x,u n ) zeB(0,u„) 
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in fact, since the graph Q is transitive, there exits an automorphism a x , such that a x (x) = 
(0 is a fixed vertex in S). Equality (fTOl) gives 

J2 Cov(Y x ,Y y )= Cov(Y ,Y ax(y) ). 

y£B(x,u n ) y£B(x,u n ) 

Now, Lemma 1.3.2 in Godsil and Royle [?] yields that d(x, y) = d(a x (x), a x (y)) = d(0, a x (y)). 
From this we deduce that y E B(x,u n ) if and only if a x (y) G B(0,u n ). From above, we 
conclude that, 

Cov{Y x ,Y y )= Gov(Y ,Y ax(y) ) = ^ Cov(Y ,Y z ), 

y&B(x,u n ) a x (y)GB(0,u n ) zeB(0,u n ) 

which proves ( |53l l. Consequently, 

' z£B(0,u n ) 

The last equality together with the first limit in (Pf6l) and (H71) . ensures 

lim T 2>n = VCov(y ,n). (54) 

n— >+oo ' » 

The second conclusion of Proposition 17.21 is proved by collecting the limits (1501 ). ( l52i l and 

(EU). □ 
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